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In the context of the dynamical mean-field theory of the Hubbard model, we identify microscopically 
an order parameter for the finite temperature Mott endpoint. We derive a Landau functional of 
the order parameter. We then use the order parameter theory to elucidate the singular behavior of 
various physical quantities which are experimentally accessible. 



PACS Numbers: 71.30.+h, 71.10.Fd, 71.27. +a 

When the strength of the electron-electron interaction 
U is increased compared to the bare bandwidth 2D, a 
metal-insulator transition (MIT) occurs §. This phe- 
nomenon, known as the Mott transition, can take place in 
the absence of magnetic long-range order, and is still an 
outstanding problem in condensed-matter physics. From 
a theoretical point of view, a difficulty is the absence 
of an obvious order parameter to systematize the criti- 
cal behavior of the observable quantities when the metal 
insulator transition is not accompanied by the onset of 
magnetic long range order. These issues are experimen- 
tally relevant to systems such as V2O3 and Ni(Se,S)2 and 
are the subject of intensive experimental study. |^| 

In recent years, great progress has been made by us- 
ing the dynamical mean-field theory (DMFT) Q. This 
framework describes both paramagnetic metallic and 
paramagnetic insulating phases. The U-T phase diagram 
(T is the temperature) of the frustrated Hubbard model 
in the limit of large lattice coordination is qualitatively 
similar to that of the V2O3 and Ni(Se,S)2 systems: A 
first-order phase-transition line ends in a second-order 
critical point, henceforth referred to as the Mott critical 
point, which is the main focus of this letter. We will 
use this framework to address the fundamental questions 
raised in the previous paragraph. 

There are two earlier qualitative ideas as to what 
should be the order parameter to describe the physics 
around the finite temperature Mott point. One idea is 
to connect the order parameter to the notions of " metal- 
licity" or coherence. It can be traced back to the early 
paper of Brinkman and Rice [Q and is captured in a slave 
boson formalism where the metallic state has a non zero 
expectation value of a Bose field which describes the co- 
herent propagation of one particle excitations. |^] In a 
very different picture, Castellani et al. viewed the metal 
as a liquid rich in doubly occupied sites, and the insulator 
as a liquid with few doubly occupied sites. The metal to 
insulator transition is viewed as a condensation of dou- 
bly occupied sites, and the order parameter is related 
to the Blume- Emery- Griffith model S. The Landau ap- 
proach presented here provides a synthesis of these ideas. 
It bridges naturally between a picture based on one par- 



ticle excitations and a picture based on local collective 
excitations (or double occupancies). In agreement with 
Castellani et al. we find that the Mott transition has 
indeed an Ising-like character. On the other hand, we 
obtain a complementary description in terms of the one 
particle spectral function reminiscent of the slave boson 
picture. A simple and clear description of the critical 
behavior near the critical point emerges. It allows us to 
systematically derive the critical behavior of any observ- 
able quantity and to relate its non analytic dependence 
on T and U to that of the order parameter. Our results 
should be also of help in resolving some controversies on 
the solution of the Hubbard model in infinite dimensions 
@i[§ by providing a theoretical framework in which to 
analyze numerical results on the finite temperature Mott 
transition. It can also be used to analyze results of pho- 
toemission and optical conductivity experiments. 

For simplicity, we focus on the single-band Hubbard 
model at half-filling, 

(y> » 

The first term describes the hopping between nearest 
neighbors on a lattice with coordination number z. The 
corresponding half bandwidth is our unit of energy, 
D = 2t = 1. The second term is an on-site inter- 
action suppressing double occupancies by imposing an 
energy cost U on each one. In the limit of infinite di- 
mensions, z — > 00, this model can be mapped onto a 
single-impurity Anderson model (SIAM) supplemented 
by a self-consistency condition. We adopt a semicircular 
density of states, which is realized on the Bethe lattice. 
The dynamical mean-field equations can be obtained by 
differentiating the Landau functional 

F LG [A] = -T£^!^ + F imp [A], (2) 

n 

with respect to the hybridization function A(iuj n ) of 
the SIAM, which has the meaning of a Weiss field. 
ibj n are fermionic Matsubara frequencies, while i r i mp [A] 
is the free energy of the SIAM, given by the action 
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Simp = SiocIA = 0} + J2a,n ft {iu n ) A(iw n ) f a {iu n ) . Here, 
5i oc [A = 0] is the action of the local / level with the 
hybridization set to zero. The first term in Eq. (|^) 
is the cost of forming the Weiss field A(ioj n ) around 
a given site, while the second one is the free energy 
of an electron at this site in the presence of the Weiss 
field. Using the Green's function of the SIAM, G(ioj n ) — 
(l/2T)SFi mp /SA(iuj n ), the mean-field equation reads 



t 2 SF hG [A] 
2T 8A(iuj n ) 



t 2 G(ico n )[A,a}- A(iuj n ) =0. (3) 



Here, a = (U, T) comprises the control parameters. This 
Landau approach was used to describe the energetics 
of the Mott transition at zero temperature ||. We 
will show that near the finite temperature Mott point, 
the Weiss field has a singular dependence which can be 
parametrized by a single number which assumes the role 
of an effective order parameter for this transition. 

As in Landau theory, we assume that a finite temper- 
ature transition exists, and derive a complete descrip- 
tion of the critical behavior near the transition as follows: 
First, we expand the mean- field equation ([}]) around the 
critical point, a c — (U C ,T C ), up to third order in the 
deviation of the hybridization function from its value at 
the critical point, 8 A — A{a c + 8a) — A(a c ), and to 
first order in 8a — (U — U C ,T — T c ). This expansion is 
well-behaved because the impurity model at finite tem- 
peratures depends smoothly on a and 8A(iuj n ). In order 
to carry out this expansion it is convenient to define a 
fluctuation matrix 



M„ 



S 2 F LG [A] 



2T SA(iLu n )8A(iuj m ) 



(4) 



critical point 



M nm has the form — 8 nm + K nm , where K nm is the 
Fourier transform of a kernel K(t, t') which is propor- 
tional to the connected correlation function of an oper- 
ator 0(r) = ffduf + (u + T)f(u), < 0(t)0{t') > - < 
0(t) >< O(r') > where the average <> is calculated 
with the action of an Anderson impurity model. It is 
well known that the correlation functions of the Ander- 
son impurity model are bounded, and therefore the Ker- 
nel K is square integrable J ^ J^drdr'lK^, t')\ 2 < oo. 
Therefore it K nm is a Fredholm operator which and has 
a discrete spectrum of eigenvalues which we labeled by 
the index I. Q 

At half-filling, particle-hole symmetry guarantees that 
the order parameter A(iw) is odd and wholly imaginary. 
Accordingly, the fluctuation matrix is real and symmet- 
ric and has real eigenvalues mi belonging to eigenvectors 
4>i{iuj n ) which can be chosen to be purely imaginary and 
to form an orthonormal basis. The critical point, in this 
description of the problem, is signaled by the appearance 
of a single zero eigenvalue, m = 0, which indicates the 
occurrence of a simple bifurcation. 



Next, we represent 8 A in the eigenbasis of the matrix 
(|f), 8A(iu n ) = J^iVi't'iii^n), where all r\i are real. By 
projecting the mean-field equation (|^) onto the eigenbasis 
4>i , we obtain an equation of the form 



F, 



( %o 3 = 



(5) 



which holds for all /. f[°' is of order 8a. i*j and F^ 

have Taylor expansions in the r]j=£o, where f[ starts 
with the linear order. We solve Eq. (0) iteratively for all 
rji^a to obtain rji^o = ai + cii] 2 + dirf^. Here, a; is of first 
order in 8a, (which assures us that the leading singular 
dependence of the spectral function is proportional to 0o) 
further corrections have the form bi t/q with bi also of or- 
der 8a. By inserting this expression into the I = case of 
Eq. (||), we derive an effective equation for the zero-mode 
amplitude r/Q. We can think of 770 as the soft mode near 
the transition and the r)i=£o as massive modes. The elimi- 
nation of the massive modes renormalizes the coefficients 
of the effective action for the soft mode. In the resulting 
cubic equation for 770, we eliminate the quadratic term by 
shifting i]q by an appropriately chosen linear function in 
8a, 7] = r/ +consti x (T — T c ) +const 2 x (U — U c ). Close 
to the critical point, 77 and rjo are dominated by non an- 
alytic terms and are therefore essentially equal. We thus 
obtain an equation of state without quadratic term in ty. 



pn + erf = h. 



(G) 



Here, all quantities are real. 

As in Landau theory, a microscopic calculation of the 
Landau coefficients (p,c,h) is difficult. However we can 
extract exact information about the critical behavior 
from the knowledge that they are smooth functions of 
the control parameters, i.e. c is finite at the critical 
point, whereas p and h are linear functions of 8a, h — 
h 1 (U-U c ) + h 2 (T-T c ) a,ndp = Pl {U-U c )+p 2 (T-T c ). 
As a consequence, r\ has a singular dependence on U and 
T near the critical point. At U — U c , and for T near T c , 

7 1 {U c ,T)^sign{h 2 /c)sigii(T-T c )\T-T c \ 1 / 3 . (7) 

The mean-field equation (||) describes the Mott tran- 
sition close to the critical point in terms of the order 
parameter rj. In this form, the analogy with the liq- 
uid gas transition is evident. The Mott transition takes 
place on the line in the U-T plane where h vanishes and 
the system has full Ising symmetry. The critical point, 
(U c , Tc), divides this line into two half-lines. On the half- 
line where T < T c , there are two solutions, rj = ±y[p/c|. 
We will see later that r\ parametrizes the strength of the 
quasiparticle resonance of the single-particle spectrum 
(see Fig. ||) . A positive or negative " field" h increases or 
decreases this component of the spectral function, respec- 
tively. The field h decreases when U or T is increased, 
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because either increase eliminates the metallic coherence 
and thus reduces the value of rj. We have used the sign 
convention whereby is positive. 

We now turn to various consequences of our construc- 
tion. From Eq. (^]), we can obtain the shape of the coexis- 
tence region near the critical point, where two solutions of 
the mean field equations coexist. It is centered symmetri- 
cally about the h = line, and its width along T = const 
lines, AC/, scales with (T c — T) 3 / 2 . The constant of pro- 
portionality is given by |/ii|)[(p2 — pih,2/hi)/3] 3 / 2 . 

An important quantity which is measured in numer- 
ical simulations is the double occupancy. It is con- 
nected to our order parameter rj as follows: (d) — 
(T/U) E„{[(iw n + n)G{iw n ) - l]e iw « 0+ - t 2 G(iuj n ) 2 } = 
(d) c + c[ d \] + c^r/ 2 . In this expansion about the criti- 
cal point, we have only retained the leading and next to 
leading nonanalytic terms responsible for the critical be- 
havior. The susceptibility \ = d(d)/dU diverges at the 
critical point. For example: 

X (U, T c ) ~ (4 rf) /3)sign(/i 1 /c)|/ ll / C | 1 / 3 |C/ - U c \- 2 '\ (8) 

The double occupancy is related to the magnetization 
by the identity ((rif — n\) 2 ) = 1 — 2(d). The magnetic 
response will therefore also exhibit nonanalytic depen- 
dences on the control parameters. 

There has been several numerical studies of the finite 
temperature Mott transition in this model. The Lan- 
dau approach predicts the functional dependence of var- 
ious quantities near the transition, and therefore the ex- 
pressions derived in this paper, are useful for interpret- 
ing the numerical work. To illustrate how our approach 
sheds new light on previously obtained numerical data 
we compare in Fig. |l| the results for the double occu- 
pancy (d) obtained within the IPT and QMC calcula- 
tions with At = 0.5/ D, after carrying out the shifts and 
the rescaling described in the figure caption. Within the 
statistical errors of the QMC calculation, the agreement 
is excellent. This surprising result is consistent with the 
Landau theory: different approximations for the solution 
of the impurity model reduce to the same Landau theory 
near the critical point, but with different values of the 
Landau coefficients. Therefore, with a suitable rescaling, 
the results near the critical point should agree with each 
other, and with a fit based on the Landau theory which 
is shown in the red line in figure 1. 

Small changes in the values of At result in shifts of U c , 
T c , and (d) at criticality, but does not change the form 
of the critical behavior. We also note that the critical 
slowing down which has been observed in the iterative 
solutions of the mean field equations are a direct conse- 
quence of the presence of the soft mode r\ described in 
the Landau approach. 

From our construction it is clear that 77 provides the 
leading non analytic behavior of the Weiss field. In or- 
der to get a better feeling for its physical significance 
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FIG. 1. Double occupation (d) as a function of U for 
different temperatures. The thin black lines denote IPT 
results for T IPT = 0.0469,0.05,0.052,0.056 (top to bot- 
tom). The thick red lines are a fit to the IPT data us- 
ing the LG theory. The circles are QMC data obtained at 
Tqmc = 1/40, 1/35, 1/32, 1/25 § . The IPT results where 
shifted by a constant -0.07 along the U axis and by -0.003 
along the (d) axis. The curves for the 3 larger temperatures 
are above T c and the lowest temperature ones (2 branches) 
are just below. The inset shows the scaling of the reduced 
temperatures [T — T c ]qmc versus [T — TJipt- 



we have to understand how it can be probed experimen- 
tally. Since the order parameter is closely related to the 
amplitude of the quasiparticle peak, photocmission is an 
ideal tool to probe the temperature and pressure depen- 
dence of the order parameter near the critical point. This 
experimental technique, in the angle integrated mode, 
would also measure the convolution of the Fermi func- 
tion with the analytically continued eigenfunction of the 
zero mode, Im(f>o(iuj n = uj — iS). To visualize the shape 
of the spectral function near the critical point we must 
resort to calculations based on analytic methods such 
as IPT. The inset of Fig. [2] shows the spectral function 
very near the critical point, computed within the IPT. 
It illustrates how the compromise between metallic and 
insulating features is realized. A finite 77, depending on 
its sign, adds or subtracts spectral weight to the coherent 
low energy feature immersed in a constant backround in 
between the Hubbard bands. The zero mode is seen to 
affect mainly the low-energy part of the spectrum, which 
determines whether the system is metallic or insulating. 
The strong temperature dependence has been noticed in 
previous theoretical and experimental studies. (!(]] Its ori- 
gin and connection to an order-parameter description of 
the Mott transition, however, had not been recognized 
until now. In the main panel of Fig. || we display the 
height of the quasiparticle peak Aq = iA(i0 + )/Trt 2 , for 
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the finite-temperature Mott point. The leading non an- 
alytic behavior of other physical quantities can be ob- 
tained along similar lines, i.e. by recognizing their cou- 
pling to the order parameter. This involves a few coeffi- 
cients, (i.e. the c^'s) which depend on the observable 
(and on the approximation method) and, as in Landau 
theory, should be taken as parameters. The dependence 
on temperature and on pressure is completely determined 
from the temperature or pressure dependence of the order 
parameter that follows from Eq. (g). ACKNOWLEDG- 
MENT This work was supported by NSF 95-29138. E.L. 
was partially supported by the Deutsche Forschungsge- 
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FIG. 2. The density of states at the Fermi energy 
p(0) = Ao as a function of temperature in the critical re- 
gion (U = 2.46316 « U c )- The singular behavior of the slope 
at T c ~ 0.046897 can be clearly appreciated. The inset shows 
the variation of the spectral function for U — U c in the vecin- 
ity of T c : dashed line for T - T c /T c = -0.00025, solid line for 
T-T c /T c = 0.00006 and dotted line for T-T c /T c = 0.00049. 



U ~ U c , as a function of temperature in the vicinity of 
T c . The rapid variation seen in the figure is consistent 
with the form Aq = Aq c + 77 + c^rj 2 with coefficients 

(A) 

c\ independent of U and temperature. 

Optical techniques are probably the best tool available 
to test the predictions of our theory. For instance, one 
may consider the integral of the optical conductivity up 
to some cuttoff, N c g(T). Since the optical conductivity 
in infinite dimensions is directly expressed in terms of 
the single-particle Green's function, N c g(T) must also 
exhibit the singular temperature dependence near the 
transition. We would therefore expect the temperature 
variation of this quantity to be most visible for a rel- 
atively small cuttoff, displaying a rapid variation with 
T similarly as for Aq. Since the singular dependence 
arises from the order parameter 77, it should be possi- 
ble to fit the Drude weight by N eS (T) = N cS (T c ) + 
c[ N) r](T) + c { 2 N) r] 2 (T). N cS (T) has recently been mea- 
sured in NiS2- x Se K [[nj, the observed strong temperature 
dependence of the effective number of carriers is consis- 
tent with our predictions. 

In summary, we derived an order parameter descrip- 
tion of the Mott transition near its critical point in the 
U-T plane. We showed that the critical behavior in prox- 
imity to this point is governed by an Ising-like Landau 
functional and is present in a large number of observ- 
able quantities. We predict that any physical quantity 
which is sensitive to the single-particle spectrum exhibits 
singular dependences on the control parameters close to 
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